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ABSTRACT 

Patterns  are  defined  to  be  stable  stationary  nonconstant  solutions 
of  the  equations  of  reaction  and  diffusion.  Several  approaches  are  used 
to  show  the  existence  (or  nonexistence)  of  patterns  depending  on  one 
variable  and  defined  on  the  entire  real  line.  For  a scalar  equation,  it 
is  shown  that  there  are  essentially  no  patterns.  For  a system,  small 
amplitude  patterns,  larger  amplitude  "peaks",  and  larger  amplitude 
"plateaus"  are  treated.  In  all  cases,  stability  is  an  important  considera- 
tion. Applications  to  ecology  and  biophysics  are  mentioned. 
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I.  Introduction  • rt 
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( u and  f are  n-vectors,  D a matrix,  often  taken  to  be  diagonal 
and  nonnegative)  have  been  the  object  of  a considerable  number  of  studies 
in  recent  years,  principally  because  of  their  actual  and  potential 
applicability  to  a variety  of  problems  in  population  dynamics,  biophysics, 
chemical  physics,  and  chemical  engineering. 

Here  we  concentrate  on  the  question  of  existence  of  stationary 
"patterns",  which  we  define  to  be  stable,  stationary,  nonconstant, 
bounded  solutions  of  (1).  (We  do  not  consider  "moving  patterns"  here, 
such  as  those  studied  by  Howard  and  Kopeil,  Winfree,  and  others.) 

Interest  in  such  stationary  solutions,  sometimes  called  dissipative 
structures,  was  aroused  by  the  work  of  Turing  [ 27  | in  the  19  50's,  and  by 
that  of  Gmitro  and  Scriven  [ 1 5 ] , Prigogine  and  Nicolis  [ 22  ) and  others 
in  the  1960's  and  later.  Particular  interest  in  them  has  been  occasioned 
by  their  possible  role  in  reflecting  the  corresponding  phenomena  of 
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The  equations  of  reaction  and  diffusion 
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pattern  iormation  in  developing  organisms  and  in  ecological  communities. 


However,  it  should  be  emphasized  that  for  most  contexts  in  which  (1) 
has  been  used  to  model  a biological  phenomenon  (signal  dynamics  on 
nerve  axons  excepted),  the  model  is  highly  idealized  and  accounts 
imperfectly  for  only  part  of  the  processes  actually  occurring  to  produce 
the  phenomenon. 

Except  for  remarks  at  the  end  of  the  paper,  we  restrict  attention 
to  solutions  depending  on  only  one  space  variable.  This  is  necessary 
for  some  of  our  results,  merely  convenient  for  others.  Possible  exten- 
sions to  many-dimensional  problems  are  noted  in  section  IV. 

We  are  particularly  interested  in  patterns  which  are  not  elicited 
by  boundary  effects,  and  for  this  reason  further  focus  attention  on 
patterns  defined  on  the  entire  real  line.  Again  in  the  last  section,  some 
discussion  will  be  given  about  corresponding  results  on  bounded  domains. 

The  above  definition  of  pattern  requires  that  it  be  stable,  with 

reference  to  the  evolution  system  (l).  For  some  results  given  here, 

this  means  C°-stability,  which  we  define  as  follows:  a bounded 

stationary  solution  <p(x)  of  (1)  is  C°-stable  if,  given  any  e > 0,  there 

is  a 6 such  that  every  solution  u(x,t)  of  (1)  defined  for  x t R,  t > 0, 

satisfying  I u(  • , 0)  - ^ | = sup  |u(x,0)  - <p(x)  I < 6 also  satisfies 

x t R 

I u(  • , t ) - <p  1 0 < e for  all  t > 0 . 
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In  other  cases,  we  are  unable  to  prove  C -stability,  and  instead  give 
stability  arguments  based  on  an  analysis  of  the  spectrum  of  the  linearization 
of  the  right  side  of  (l)  about  the  stationary  solution  in  guestion.  In  still 
other  cases,  we  give  less  rigorous  heuristic  arguments  favoring  stability. 

Variations  of  the  above  reaction-diffusion  problems,  such  as  those 
obtained  by  allowing  the  diffusion  terms  to  be  nonlinear,  and/or  first  deriva- 
tive terms  to  enter  the  equation,  are  important  in  some  applications,  but  will 
not  be  discussed  here.  Othmer  [ 21  ) has  discussed  (among  other  things) 
standing  oscillation  solutions  of  (l),  which  bear  the  same  relation  to  x- 
independent  oscillatory  solutions  as  do  ours  to  constant  so.utions. 

The  paper  is  concerned  with  the  mathematical  analysis  of  patterns; 
experimental  and  numerical  simulation  results  are  not  mentioned  except  in 
passing.  Nor  do  we  go  into  much  detail  about  applications  in  the  various 
fields  mentioned  above.  For  an  extended  discussion  of  the  relevance  of 
patterns  in  ecology,  see  (for  example)  the  article  by  Levin  [ 18]  . 

The  plan  of  the  paper  is  as  follows. 

I.  Introduction 

II.  The  scalar  case,  n = 1 

III.  Systems,  n > 1 

A.  Small  amplitude  patterns 

B.  Larger  amplitude  patterns 

1 . Peaks 

2.  Plateaus 

IV.  Discussion 
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Throughout,  we  assume  that  f is  a continuously  dillerentiaule 
function  of  its  argument. 

I gratefully  acknowledge  helpful  discussions  with  D.  H.  Sattinger 
and  with  K.  Kirchgassner,  who  convinced  me  that  there  exist  quasiperiodic 
solutions  of  (14),  and  that  only  the  cases  shown  in  Fig.  4 occur. 
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II.  The  scalar  case,  n 1 

In  this  case,  D is  a scalar,  which  we  take  to  be  unity.  The 
equation  of  one-dimensional  patterns,  then,  is  simply 

(2)  u + f(u)  = 0,  x t R . 

xx 

This  equation  is  well-known  and  easily  solved  by  a phase-plane  analysis 
or  by  the  argument  below.  Furthermore,  the  stability  of  its  bounded 
solutions  can  be  tested  by  means  of  a technique  due  to  Aronson  and 
Weinberger  [ 1 ) . 

The  equation  has  a first  integral,  obtained  by  multiplying  by  the 
integrating  factor  u^: 

1 2 

(3)  ~ (u  ) + F(u)  = E = const  , 

L.  X 

where 

u 

F(u)  = f f(s)ds  . 

0 

If  one  draws  the  graph  of  the  "potential"  V = F(u),  as  illustrated 
in  Fig.  1,  then  there  is  a one-one  correspondence  between  the  nonconstant 
solutions  of  (2)  (modulo  translation  of  the  independent  variable  and 
reversal  of  its  sign,  which  are  always  possible)  and  horizontal  line 
segments  in  the  (u,V)  plane  whose  finite  endpoints  lie  on  the  given 
curve,  and  which  otherwise  lie  strictly  above  the  curve.  The  projection 
of  such  a segment  onto  the  u-axis  is  the  range  of  the  corresponding 
solution,  and  the  ordinate  of  the  segment  is  the  constant  E (its  "energy") 
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Fig.  1.  The  potential  function  V = F(u).  Horizontal  lines 
(and  points)  such  as  a - f correspond  to  solutions  of  (2). 


in  (3).  Given  such  a segment,  the  solution  u(x)  can  be  obtained  by 
further  integrating  (3).  On  each  interval  of  monotonicity,  we  obtain 

i u i 

(4)  x - xQ  = ± 2~ i f (E  - F(s)f  *ds  , 

U0 

where  uQ  is  an  arbitrary  number  in  the  interior  ot  the  (known)  range, 
and  xQ  is  arbitrary.  This  excludes  the  constant  solutions  of  (2); 
but  they,  of  course,  can  be  identified  with  the  points  of  zero  slope  in 
Fig.  1,  since  then  F'(u)  : l(u)  - 0. 


Fig.  1 illustrates  six  solutions:  (a)  and  (d)  are  constant  ones; 

(b)  is  a periodic  solution;  (c)  is  a solution  attaining  a minimuiu  ut  a 

single  finite  value  of  x,  and  approaching  its  supremum  as  x - ±ae; 

(e)  is  a monotone  solution  approaching  different  limits  as  x - ±»; 

and  (f)  is  a solution  bounded  from  below,  but  unbounded  above.  These 

properties  can  be  deduced  from  (4).  For  example,  in  cases  (c)  and  (e), 

the  tact  that  lx  I as  u approaches  its  supremum  u follows 

m 

2 

from  the  fact  that  E - F(s)  < C(u  - s)  for  s near  u , so  that 

— m m 

the  integral  in  (4)  is  unbounded  as  u f u . In  the  case  (b),  a similar 
analysis  shows  that  the  maximum  and  minimum  are  obtained  for  finite 
values  of  x.  Then  the  uniqueness  of  initial  value  problems  for  (2), 
together  with  the  sign  reversibility  of  x,  show  the  solution  must  be 
periodic. 

With  all  possible  solutions  of  (2)  known,  we  now  consider  their 
C°-stability  as  solutions  of  (1),  which  in  our  present  case  is 
(S)  u - u - f(u)  = 0 . 

t XX 

The  result  is  that  very  few  are  stable. 

Lemma  1:  Every  bounded  nonconstant  solution  of  (2)  which  attains  a 
maximum  or  minimum  at  a finite  value  of  x is  C^-unstable. 

Proof:  Let  <r  be  such  a solution.  Assume  it  has  a minimum  <r  = m 
at  x = 0 (the  argument  for  the  case  of  a maximum  is  the  same).  Then 


it  corresponds  to  a segment  I,  such  as  that  in  Fig.  2,  whose  left 


endpoint  is  on  the  curve  V = F(u)  at  a point  u = m such  that 
F'(m)  < 0.  Let  E be  the  ordinate  of  the  segment  I (the  energy  of  </>), 

and  for  positive  6 let  Ic  be  the  segment  with  energy  E + 6 

0 

overlapping  I.  Let  <p  be  the  corresponding  (possibly  unbounded) 

o 


Fig.  2 


solution  of  (2).  For  small  enough  6,  it  will  attain  a minimum  rru  < m, 
where  F(mc)  = E + 6.  Furthermore 

O 

(6)  lim  mc  = m . 

610  6 

By  translation  of  the  independent  variable,  if  necessary,  we  may 

assume  that  <*>  also  attains  its  minimum  at  x = 0. 

0 

Let  X.  be  the  largest  interval  on  the  x-axis  containing  the 

O 

origin,  on  which  to  (x)  <tp(x).  Since  is  longer  than  I at  both 


-8- 


ends,  </>  will  assume  values  greater  than 
o 

is  finite.  Let 


sup  </>(x). 
x 


Therefore  X, 
0 


<Mx) 


^(x),  x t Xf)  ; 

<*>(x),  x \ X6  . 


We  shall  need  the  following  fact,  whose  proof,  being  technical, 
is  deferred  until  later. 

Lemma  2:  lim  I ip  - <p  | = 0 . 

6 10  0 

Locally,  <p  is  the  minimum  of  two  solutions,  so  will  itself  be 

a supersolution  of  (5).  Let  u(x,  t)  be  the  solution  of  (5)  satisfying 

u(x,0)  = <*>(x).  By  a result  of  Aronson  and  Weinberger  [ 1 ] , u(x,  t) 

is  decreasing  in  t,  pointwise  in  x.  (The  proof  in  [ 1 ] was  for  the 

case  v>  = const,  but  works  as  well  in  our  case. ) 

Let  [ 0,  T)  (T  < oo)  be  the  maximum  interval  of  existence  of  u. 

By  the  decreasing  property,  lim  u(x,t)  s i|>(x)  exists  and,  where  i);  > -00, 

t t T 

is  a solution  of  (2)  (see  [1]). 

If  4j  is  a solution  for  all  x,  then  it  must  satisfy  vp(x)  < sup  ip 

and  must  attain  values  < mc  < m.  This  follows  from  the  maximum 

0 

principle  and  the  fact  that  <p(x)  < ?(x).  But  it  may  be  seen  that  any 
horizontal  segment  or  point  in  the  u - V plane  corresponding  to  such 
a solution  must  lie  entirely  to  the  left  of  the  value  u = m - r^,  for 
some  positive  n independent  of  6.  Thus  for  each  x, 

(7)  u(x,  t)  < m - r, 

for  large  enough  t where  u is  defined. 
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On  the  other  hand,  if  « is  not  a solution  for  ali  x,  then 
v -ac  for  some  x,  so  in  any  case,  ( i ) holds  for  at  least  some  values 

of  x and  large  enough  t.  (It  can  be  shown,  in  fact,  that  either 
4;  £-oo,  or  4,  = const,  and  that  u - 4,  as  t - « uniformly  on 
bounded  sets.  But  we  shall  not  use  these  results.) 

Now  let  u be  any  solution  of  (5)  satisfying  u(x,0)  < <p(x).  By 
the  maximum  principle,  u(x,  t)  < u;x,  t),  so  there  are  some  values  oi 
x for  which  u(x,  t)  < m - g for  all  large  enough  t where  u is 
defined.  For  such  t,  we  have  I u(  • , t)  - <p  lQ  > g > °-  But  in  view 
of  Lemma  2,  we  may  choose  !u(-,0)  - I q arbitrarily  small,  therefore 
c is  unstable,  and  the  theorem  is  proved. 

Besides  the  unstable  ones  covered  in  the  theorem,  the  only 
nonconstant  bounded  solutions  of  (2)  are  the  monotone  ones,  correspond- 
ing to  (e)  in  Fig.  1.  These  can  be  thought  of  as  zero  velocity  travelling 
wave  solutions  in  the  sense  of  ( 1 3],  where  their  global  stability  was 
proved.  On  the  other  hand,  their  existence  depends  on  a special 
property  of  the  curve  F in  Fig.  1,  namely  that  two  adjacent  local 
maxima  have  the  same  height.  A slight  change  in  f will  destroy 
this  property,  so  the  monotone  solutions  are  structurally  unstable. 

This  exhausts  the  possible  bounded  solutions  of  (2),  except 
for  the  constant  ones.  Though  they  are  not  candidates  for  being  patterns, 
we  include  a discussion  of  their  stability  for  the  sake  of  completeness. 

As  it  turns  out,  a constant  solution  u 5 u^  of  (2)  is  C -stable  if 
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and  only  if  it  is  stable  as  a solution  of  the  equation 


1«)  77  - f(u)  . 

dt 

Suppose  first  that  uQ  is  unstable  as  a solution  of  (8).  Then 
there  exists  a neighborhood  T of  uQ,  and  solutions  u(t)  of  (8) 
with  u(0)  arbitrarily  close  to  u^,  but  with  u(t)  f 7>  for  some  t. 

Such  a function  u(t)  is  also  a solution  of  (5),  however.  Therefore  u^ 
is  C -unstable  as  a solution  of  (5)  as  well.  Next,  suppose  uQ  is 
staple  as  a solution  of  (8).  Then  for  each  e > 0,  there  is  a 6(e)  > 0 
such  that  solutions  of  (8)  with  lu(0)  - u^  I < b exist  for  all  time  and 
satisfy  I u(t)  - u I < e . Let  u(x,  t)  be  a solution  of  (5)  with 
I u(  • , 0)  - Ug  I < 6 . By  the  maximum  principle,  it  is  bounded  pointwise 
above  and  below,  respectively,  by  u (t)  and  u (t),  which  are  solutions 
of  (8)  assuming  initial  values  u±(0)  = uQ  ± 6 . Therefore  lu(*,t)-u  ! < e, 
and  u(x,  t)  is  stable. 

Altogether,  we  have  the  following  theorem. 

Theorem:  The  only  stable  solutions  of  (2)  are  the  strictly  monotone  ones, 
and  the  constant  ones  which  are  stable  as  solutions  of  (8).  Moreover, 
the  strictly  monotone  ones  are  structurally  unstable. 

Proof  of  Lemma  2:  First,  consider  the  case  when  v is  periodic.  Then 

geometrical  considerations  show  that  the  length  of  Xf  will  not  surpass 

twice  the  wave-length  of  <fi,  so  remains  bounded  independently  of  6. 

Then  the  continuity  of  solutions  of  initial  value  problems  with  respect 

to  initial  data  easily  implies  lim  lv>(x)  - v’(x)  I - lim  - «>t(x))  - 0, 

610  610 
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t 


uniformly  in  X.  . But  outside  X.,  v’ - v , so  the  conclusion  follows. 

The  only  other  case  is  when  <p  approaches  a limit  asymptotically 
as  lx!  * »,  and  v>'  > 0 for  x > 0.  We  treat  this  case  as  follows. 
By  asssumption,  r'(m)  = f(m)  < 0.  Let  k = ~ f(m)  < 0,  and  restrict  6 
to  be  so  small  that  f(u)  < k for  m,  < u < m.  Then 

x s 

*6(x)  = v6(0)  + f f ^'(4)dgds 


X s 

= *(0)  - f f f(«M|))d£ds 
o o 

1 i I 2 

> m&  + - I Klx  » 


for  x such  that  mc  < <p, 

0 ~ O 

(x)  < m.  Let 

2 

w6 

2(m  - mf )/ 1 «l , so  that 

m.  + ~ Ixluif  = m.  Then 
o 2 o 

v6(u6)  > m. 

If  we 

define  it,  as  the  first 

value  of  x > 0 at  which 

<^6(x)  = m, 

then 

tt , <u>,.  It  follows  that 

0 — 0 

(9) 

lim  tt. 

. . „ 0 

= 0 . 

6 1 0 

Let  e <p,(x  + nft).  Then  4^(0)  = m = «>(0),  and  from  (3),  for  x > 0, 

I i 1 

^ = [ 2(E  +6  - F^))]2,  whereas  = [ 2(E  - F(</>))  | 2 < [ 2(E  + 6 - F(<^)) ) - . 

It  follows  that  4^  (x)  > <p(x),  for  x > 0 . 

We  now  have,  for  x > nf  , 

<p(x)  - <Pf  (x)  ~ <P(X)  ~ <P(*  - 1T&)  + i P(X  - TT^)  - 4>&(X  - T7f  ) 

< tp(x)  - <p(x  - irf|)  . 

By  the  uniform  continuity  of  </>  and  (9),  we  have  that  <p(x)  - <p^(x)  < n I 0 
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as  6-0,  for  x >ir  , By  symmetry,  the  same  is  true  for  x < -n, 

o o 

as  well,  and  by  the  continuity  of  v5  and  <p^ , we  may  extend  the 
relation  to  the  interior  interval  | -ttc  , ttc  1 . 

o 6 

On  X we  therefore  have 

o 

l10)  l</>(x)  - <p(x)|  = <p(x)  - <pAx)  < <x  , 

whereas  outside  X^,  we  have  l^(x)  - ip(x)|  = 0.  This  completes  the 
proof. 


HI.  Systems  (n  I) 


A.  Small  amplitude  perturbations. 

Let  us  consider  systems  (1)  in  which  f(0)  0,  so  that  u s o 

is  a solution.  It  has  been  observed  in  the  past  that  this  trivial  solution 
may  be  unstable  with  respect  to  certain  perturbations  which  are  not 
constant  in  x,  while  remaining  stable  with  respect  to  constant  (uniform) 
perturbations.  Another  way  of  saying  this  is  that  the  trivial  solution 
u = 0 may  be  unstable,  whereas  the  point  u - 0 is  a stable  critical 
point  of  the  associated  kinetic  equations 

(11)  77  = f(u)  . 

dt 

In  such  cases,  the  effect  of  introducing  diffusion  and  x-dependence  is 
therefore  to  destabilize  the  system,  somewhat  contrary  to  one's  intuition 
of  diffusion  as  a stabilizing  influence.  An  argument  was  advanced  by 
Segel  and  Jackson  [ 23)  to  partially  explain  this  paradox.  It  one  can 
characterize  some  of  the  components  of  u as  "stabilizing"  in  some 
sense,  and  others  as  "destabilizing",  and  if  the  stabilizers  diffuse  more 
rapidly  than  the  destabilizers,  then  a very  small  initial  nonuniform 
concentration  of  both  (this  being  the  perturbation  of  the  trivial  solution) 
may  result  in  the  stabilizers  diffusing  away,  their  effect  thereby 
diminishing.  This  would  happen  to  some  extent  not  only  with  the  initial 
concentration  of  stabilizers,  but  also  with  stabilizers  produced  later 
through  the  reaction  process.  In  this  way,  without  as  much  counteracting 
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influence,  the  destabilizers  could  take  over,  and  the  solution  could 

grow  to  finite  size.  Tor  this  to  happen,  it  is  of  course  necessary  that 

the  diffusion  coefficients  not  all  be  the  same.  In  the  example  in  [24], 

n = 2 and  the  system  represented  predator-prey  equations,  the  predators 

being  stabilizers  and  the  prey  destabilizers.  Segel  and  Levin  [ 2S| 

have  pursued  similar  problems  and  have  shown  that  after  instability  sets 

in,  one  can  expect  solutions  of  (1)  with  small  initial  data  sometimes 

to  evolve  into  patterned  solutions. 

We  consider  the  system  (1),  and  assume  the  function  f depends 

on  a real  parameter:  f = f(u,\).  We  also  assume  that  there  is  a critical 

number  \ such  that  for  \ < X , the  trivial  solution  of  (1)  is 
cr  cr’ 

stable,  in  a sense  to  be  explained  below,  whereas  for  X > \ , it 

cr 

is  not.  The  new  instability  is  respect  to  periodic  perturbations  with 

wave  length  in  a certain  range,  bounded  above  and  below  (the  meaning 

of  this  is  given  below).  In  particular,  the  zero  solution  remains  stable 

to  constant  perturbations,  so  the  instability  is  of  the  type  indicated 

in  the  previous  paragraph.  Under  certain  additional  conditions,  the 

effect  of  this  is  that  periodic  patterns  of  small  amplitude  arise  for  X 

in  a right-hand  neighborhood  of  X .In  fact,  we  shall  show  how 

cr 

bifurcation  theory  may  be  used  to  construct  a two-parameter  family  (the 
parameters  being  amplitude  and  wave-length)  in  a neighborhood  of  \ . 

The  question  of  their  stability  will  be  taken  up  later. 


Our  assumptions  need  to  be  made  precise.  The  source  term  f 

will  be  assumed  twice  continuously  differentiable  for  u near  the 

origin;  we  write  A for  the  Jacobian  ^(0,X).  Tor  simplicity,  we 
\ Bu 

assume  the  dependence  of  on  X is  linear: 

some  matrices  A and  B.  The  linear  problem  associated  with  (1), 

V = DV  + AV  + XBV , 
t xx  ’ 

has  solutions  of  the  form  V = <1*  exp[  ikx  + <rt)  for  arbitrary  real 
wave  number  k.  Here  will  be  an  eigenvector  of  the  matrix 

H(X,  p)  = -pD  + A + XB  (p  = k^)  , 
with  corresponding  eigenvalue  cr. 


A^  = A + \ B,  for 


Assumption  1:  There  is  a number  \ and  a number  6 > 0.  such  that 
cr 

for  i\  - \ i < 6 and  for  all  p > 0,  H(\,  p)  has  a unique  eigenvalue 

with  largest  real  part,  and  it  is  real  and  simple.  Denote  it  by  Oj(X,p), 

and  let  I = { p : <x  (\ , p)  > 0 } . Then 
\ i 


0 for  \ - b < \ < \ , 

cr  cr 

1 pQ } for  X - X , where  pQ  > 0, 

a bounded  interval  of  positive  length  for  X < X < X + 6 . 

cr  cr 


This  assumption  about  o (X,p)  is  depicted  in  Fig.  3. 

For  simplicity  we  take  X ^ - 0 from  now  on  (except  in  the 
example  below). 

It  follows  from  Assumption  1 that  the  matrix  HQ  = H(0,  pQ)  has 
zero  as  a simple  eigenvalue.  Its  nullspace  is  therefore  spanned  by  some 
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Fig.  3.  Typical  dependence  of  the  largest  eigenvalue  on 

V and  p,  when  Assumption  1 is  valid.  The  curves  are 
for  \ constant. 


real  nullvector  and  the  nullspace  of  the  adjoint  matrix  is  spanned 

by  some  real  vector 
Assumption  2:  B®  4 0. 

Example:  Let  n = 2, 


** 

»— 

— 

0 

-2 

1 

0 

1 0 

A = 

2 

_ 1 
X 

, B = 

0 

0 

, and  D = 

0 2- 

Then  Assumption  1 is  satisfied  with  p0  = n/2  - ^ , V = ^(4V2  - 1). 
In  fact,  the  eigenvalues  <r  of  H(\,p)  are  the  roots  of  the  equation 
o*  + o(  3p  + 1 - \)  + Q(p,\)  = 0,  where  Q(p,  X)  (p  - \)(2p  + 1)  + 4. 
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It  can  be  shown  that  for  small  ! \ — \ I,  the  two  roots  are  real  and 

cr 

distinct  for  all  p > 0,  and  the  lesser  of  the  two  is  negative  and 
bounded  away  from  zero.  The  greater  one,  u (p, \),  can  be  shown  by 
the  quadratic  formula  to  satisfy 

Oj  ^ 0 according  as  Q(p, \)  - 0 . 

But  Q(p, \)  is  quadratic,  and  has  a minimum  with  respect  to  p for 


each 

\.  Its  minimum  is  0 (at  p - 

P0)  for 

II 

o 

is  positive 

when 

X < X , 
cr 

and  is  negative  when 

X > \ . 

cr 

Thus  the 

set  of  values 

of  p 

for  which 

o >0  is  empty  for 

\ < \ , 
cr 

and  is  a 

finite  interval 

for  \ 

> X . 

cr’ 


The  basic  result  is  the  following. 

Theorem : Assume 

(12)  f(u,\)  = (A  + \B)u  + G(u)  , 
where  G is  continuously  differentiable  near  u = 0,  and 

(13)  |G(u)|  < C|u|2 

for  iuI  < 1.  Let  A,  B,  and  D satisfy  Assumptions  1 and  2.  Let 
kQ  \'pQ  be  the  "initial"  wave  number.  Then  there  exists  a two- 
parameter  family  of  periodic  functions  u(x;c,k),  and  a scalar  \(e,  K), 
both  continuous  in  e and  « and  defined  for  small  lei,  Ik!,  such  that 
(i)  u and  X.  satisfy 

(14)  Du  + (A  + \B)u  + G(u)  = 0,  x t R , 
for  each  fixed  e and  k, 
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(ii)  u is  periodic  in  x with  wave-length  (period)  Zn/k 
where  k = k t k, 

£n/k 

(iii)  j u(x;e,  k)  • cos  kxdx  = e , 

0 

(iv)  \(0,  0)  - 0,  u(x;  0,  0)  0. 

We  shall  sketch  the  proof.  The  first  step  will  be  to  give  the 
problem  a more  convenient  formation.  By  using  the  new  variable  x x/l  , 
we  normalize  the  wave-length  of  our  solutions  to  be  Ztt,  but  at  the 
same  time  introduce  an  extra  parameter  k^  = p into  the  equations  (14). 
which  become 

(15)  pDu--  + (A  + \B)u  + G(u)  = 0 . 

XX 

Henceforth  we  shall  omit  the  bars  on  the  x's.  Tor  p = p , X - 0,  tne 
corresponding  linear  equation 

p DU  + AU  = 0 
0 xx 

has  exactly  two  linearly  independent  solutions  bounded  for  all  x : <r  - cos  x 
and  <*>'  = -<*>  sm  x.  This  follows  from  Assumption  1.  We  normalize  our 

2it  l z 

solutions  so  that  (<r,<r)  = f Iv!  dx  ■ ■ 4 I it  1.  Similarly,  the 

0 

linear  adjoint  problem  has  two  linearly  independent  solutions  y 4-  cos  x 
and  4j'  = -*  sin  x. 

Let  Pj,  P^  be  S ^-orthogonal  projections  onto  span  {jj}  and 

span  { ip* } respectively;  for  example,  P.f  = (f,  4)4-  Also  let  Q = id  - P - P . 

1 I <- 
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Let  Z be  the  Banach  space  of  continuous  ZiT-periodic  functions  with 
values  in  Rn,  the  space  of  Revalued  functions  with  bounded 

continuous  derivatives  of  order  k,  X = Z (1  and  Y the  Z 
orthocomplement  of  <p  and  <f'  in  X.  We  norm  X and  Y by  the 
usual  norm,  and  Z by  the  C°  norm. 

Then  it  is  standard  that  the  differential  operator 

,2 

L-  P0D—  + A, 
dx 

restricted  to  Y,  is  one-one  from  Y onto  QZ,  so  that  L * exists 
as  a bounded  operator  from  QZ  to  Y. 

Since  X = span  {<?,<?' } © Y,  any  solution  in  X can  be  written  as 
u - a<p  + fi<p'  + v = 4>(a  cos  x + p sin  x)  + v,  v t Y.  Furthermore  since  (15) 
is  autonomous,  any  solution  remains  a solution  when  subjected  to  an 
arbitrary  translation  of  the  independent  variable.  In  fact,  such  a transla- 
tion, does  not  remove  v from  Y.  We  suppose  our  solutions  have  all 
been  so  translated  in  such  a way  that  p = 0;  we  then  let  e = a and 
call  it  the  "amplitude"  of  the  solution.  In  short,  we  seek  solutions 
in  the  form 

(16)  u = z{<p  + w),  w c Y . 

With  the  two  final  definitions  F(u,e)  = e *G(eu)  (recall  (13)) 
and  q = p - pQ,  we  are  ready  for  the  reformulation.  Setting  (16)  into 
(16)  and  using  the  above  definitions,  we  see  that  the  problem  is  to  find 
a function  w t Y satisfying 
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(IT)  Lw  + \B(«>  + w)  4 qD(<p  4 w)"  4 cT(<p  4 w,  e)  = 0 . 

We  rewrite  this  os  Lw  4 R(w;X,e,q)  = 0,  where  for  each  (X,e,q),  R 
is  an  operator  from  Y into  Z.  Applying  projections  Q,  P and  P 
we  see  that  (17)  is  equivalent  to  the  three  equations 


(18a) 

Lw  4 QR( w; X,  e,  q)  = 0 

(18b) 

Pj  R(w; X,  e,  q)  = 0 

(18c) 

P^R(w;X,  e,  q)  = 0 

With  the  aid  of  Assumption  Z and  the  implicit  function  theorem, 
the  first  two  equations  (18a,  b)  may  be  solved  uniquely  for  w(e,  q)  and 
Ve,q)  for  small  (e,q).  It  remains  to  be  shown  that  this  solution  will 
automatically  satisfy  (18c).  Here  we  use  the  following  simple  invariance 
property  of  our  problem.  The  operators  L,  R,  Q,  and  L *Q  preserve 
the  class  of  even  functions.  This  fact  is  immediate  for  L,  R,  and  Q, 
and  not  too  difficult  for  L ^Q.  It  follows  that  (18a-b)  may  also  be  solved 
in  the  class  of  even  functions,  so  by  uniqueness,  the  original  solution 
must  be  even.  But  P . being  orthogonal  projection  onto  a space  of 
odd  functions,  annihilates  even  functions.  Therefore  (18c)  is  satisfied, 
completing  the  proof  of  the  tneorem. 

Equations  (18)  lend  themselves  to  the  calculation  of  explicit 
approximations  to  w and  X.  as  power  series  in  e and  q.  Upon 
calculating  these  power  series  to  terms  of  order  two,  one  finds 

(19 ) X = Yje^  + Y + + higher  order  terms  , 
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where  the  coefficients  y can  be  written  explicitly  in  terms  of  various 

2 

combinations  of  the  operators  B,  D r,  L Q,  P , F(-,0),  and  first 

dx  1 

order  differentials  of  F,  acting  on  the  function  <p.  The  fact  that  there 
are  no  first  order  terms  in  the  expression  for  \ is  a consequence  o(  the 
facts  that  F(  • , 0)  is  quadratic  or  zero,  that  <p  is  a vector  multiple 

9,1 1 

of  cos  x,  and  that  _ — (0,  p ) - 0 (in  terms  of  Fig.  3,  <r  (X.  , p) 

8p  o 1 cr 

has  a maximum  at  p = p^). 

8o 

If  we  further  assume  that  r-(0,  p ) < 0 and  — ~ (0,  p ) > 0 

c.  0 (j \ u 

9p 

(which  inequalities  are  reasonable,  from  Assumption  1),  then  it  follows 
that  \ ^ > 0 . In  this  case,  it  can  be  seen  that  the  "bifurcation  diagram" 
is  generically  of  one  or  the  other  of  the  two  types  shown  in  Fig.  4.  In 
that  figure,  the  shaded  regions  consist  of  those  pairs  (X.,  e)  for  which 
a periodic  solution  u = z(<p  + w)  exists.  These  are  local  diagrams: 
e,  q,  and  X.  are  all  assumed  small. 

In  the  context  of  problem  (14),  we  have  shown  that  when  increasing 
the  value  of  X.  causes  the  trivial  solution  to  lose  its  stability,  it 
also  causes  many  new  periodic  solutions  to  appear.  There  are  two 
important  questions  now  to  be  asked:  (i)  are  there  any  other  small  bounded 
solutions  besides  the  periodic  ones  we  have  constructed:  and  (ii)  which 
c f our  new  solutions,  if  any,  are  stable?  These  are  difficult  questions, 
and  no  complete  answer  to  them  exists  at  this  time.  However,  some 
comments  can  be  made. 
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>-  x 

(a) 

Tig.  4.  The  shaded  areas  are  the  totality  of  small  (\,e)  for 
which  a periodic  solution  exists,  with  wave  number  near  k^: 

(a)  the  case  v.  = V.  “ y?(4vJ  1 > 0;  lb)  the  case  y < 0. 

4 1 £ 3 4 

Regarding  the  first  question,  we  have  constructed  all  possible 
periodic  solutions  with  wave  number  near  the  one  (kQ)  with  respect 
to  which  instability  first  sets  in.  This  is  clear  by  the  implicit  function 
theorem's  uniqueness  assertion.  However,  there  undoubtedly  exist 
quasiperiodic  solutions  as  well,  as  Kirchgassner  has  argued  in  ms  treat- 
ment of  another  type  of  model  problem  [17].  In  fact,  this  is  clearly  true 
of  (17)  in  the  case  e = 0 (the  limit  equation  for  small  amplitude  solu- 
tions). This  equation  is  linear  with  constant  coefficients,  and  tor  each 
small  \ > 0,  there  are  exactly  two  values  of  q for  wnich  it  has  even 
periodic  solutions.  In  terms  of  the  original  formulation  (14).  they 
correspond  to  two  infinitesimally  small  even  periodic  solutions,  with 
distinct  periods.  Such  solutions  of  the  linear  problem  may  be  superposed; 
linear  combinations  of  them  yield  new  quasiperiodic  solutions  (which  are 
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actually  also  periodic  with  possibly  large  period,  if  the  two  basic  periods 
are  rationally  dependent).  Each  of  the  basic  solutions  may  also  be 
subjected  to  arbitrary  translations  in  x,  so  they  are  no  longer  even. 

In  this  way,  two  more  parameters  are  introduced,  and  in  all  we  have  a 
four-parameter  family  of  quasiperiodic  solutions  of  the  linear  problem. 

For  fixed  V > 0,  one  expects  quasiperiodic  solutions  to  exist  also  for 
e > 0,  and  there  are  ways  to  construct  formal  approximations  to  them 
in  some  cases.  If  it  could  be  shown  that  a four  parameter  family  of 
them  exists,  for  some  fixed  A.  > 0,  then  it  would  follow  from  a result 
on  manifolds  of  neutral  stability  [8]  that  in  fact  there  can  be  no  other 
small  bounded  solutions  besides  those,  and  the  answer  to  (i)  would  be 
more  nearly  complete. 

As  for  the  second  question,  concerning  stability,  no  completely 
rigorous  results  will  probably  be  obtained  for  some  time.  Nevertheless, 
heuristic  results  on  linearized  stability  may  be  within  reach.  Specifically, 
on  linearizing  the  operator  on  the  left  of  (17)  about  a given  periodic 
solution  u = c(<p  + w ),  one  obtains  an  operator  with  periodic 

e,  q e, q 

coefficients,  depending  also  on  the  parameters  e and  q.  One  then 
would  look  for  all  eigenvalues  of  this  operator,  acting  on  functions, 

and  ask  for  a relation  between  e and  q which  will  guarantee  that  no 
such  eigenvalues  lie  in  the  right  half-plane  (0  is  always  an  eigenvalue). 
Such  pairs  (e,q)  would  then  yield  periodic  solutions  which  may  be 


stable  in  some  sense. 


By  a Floquet-type  argument,  it  is  sufficient  to  pose  this  eigenvalue 

icy  x 

problem  in  the  class  of  functions  of  the  form  e v(x)  with  v 2ir-periodic 
and  a real.  Then  the  set  of  eigenvalues  will  depend  on  a as  well 
as  on  e and  q,  and  one  can  restrict  attention  to  small  a,  because 
otherwise  all  eigenvalues  can  be  shown  to  have  negative  real  part. 

My  conjecture  is  that  in  case  (a)  of  Tig.  4,  only  the  solutions 
with  extremal  e are  stable  in  the  above  sense,  whereas  in  case  (b), 
none  are  stable  except  the  trivial  solution  with  \ < 0. 


III.  B.  I.  Larger  amplitude  patterns:  peaks 

There  has  been  a sizeable  amount  of  numerical  simulation  in  recent 
years  directed  toward  the  discovery  of  finite  amplitude  patterns  in  model 
reaction-diffusion  systems.  The  major  efforts  have  been  by  a group  in 
Brussels  (see,  for  example,  [ ^ | , [ 16  J ) . and  by  Gierer  and  Meinhardt  ([  14  | , 

[ 19  |,  [*10  |)  in  Tubingen.  These  efforts  have  indeed  been  successful, 
as  one  can  judge  by  reading  the  cited  papers. 

The  principal  impetus  has  been  biological;  for  example,  the  results 
bear  on  the  formation  of  morphogenetic  fields  in  developing  organisms.  The 
models  of  Gierer  and  Meinhardt  have  n - Z,  the  two  reacting  components 
being  an  "activator"  u and  an  "inhibitor"  v.  These  authors  give 
arguments  to  explain  why  patterns  may  be  expected,  if  the  inhibitor 
diffuses  more  rapidly  than  the  activator.  (There  is  some  analogy  with  Segel 
and  Jackson's  arguments  involving  stabilizers  and  destabiiizers) . 

Here  we  deal  with  the  following  "inverse"  problem:  given  a pair 
of  functions  (<r(x),  4>(x)),  find  a reaction-diffusion  system  with  n = 2, 
for  which  (</>,  i|>)  is  a pattern.  Note  that  our  theorem  in  Sec.  II 
indicates  that  the  corresponding  problem  with  n 1 does  not  have  a 
solution,  unless  </  is  monotone;  and  even  then,  the  pattern  is 
structurally  unstable.  We  prove  that  for  "single  peak"  distributions 
(<p,  ip)  of  a certain  type,  in  which  <p  and  4j  are  related  linearly,  a 
solution  of  the  inverse  problem  does  exist  (in  fact,  many  solutions 
exist).  We  also  show  that  for  "multiple  peak"  distributions,  solutions 


exist  if  the  "stability”  requirement  for  patterns  is  weakened  to  the 
statement  that  the  spectrum  of  the  linearized  operator  has  no  points  in 
the  unstable  half-plane.  It  should  be  emphasized  that  it  is  easy  to  find 
systems  for  which  a given  distribution  (</>,  4>)  is  3 solution;  but  we  must 
also  prove  that  (v,d>)  is  stable  for  the  system,  a considerably  more 
difficult  task.  The  systems  we  construct  are  of  the  activator-inhibitor 
type,  and  the  inhibitor  does  (as  in  ( 14 ) , | 19  j , ( 20  ])  di  ffuse  more  rapidly 
than  the  activator.  Our  results  are  not  likely  to  be  of  practical  importance, 
because  model  systems  (1)  which  have  been  used  in  the  past,  and 
probably  future  models  as  well,  have  other  desired  properties  besides 
merely  being  of  activator-inhibitor  type.  Thus,  the  source  functions  f 
and  g used  in  [ 14  ] , [ 19  ] , [20  J were  built  on  the  basis  of  more  or  less 
specific  morphogen  reactions.  At  the  same  time,  in  constructing  model 
reaction-diffusion  systems  with  only  two  components,  one  should  not 
attach  overriding  importance  to  having  them  mirror  specific  reaction 
networks  involving  the  two  species.  In  fact,  the  actual  mechanisms 
modelled  v/ill  involve  a large  number  of  reacting  species;  if  one  pictures 
the  reduction  to  two  as  having  been  made  through  various  pseudo-steady- 
state,  slow  reaction,  or  other  approximations,  then  the  connection 
between  the  source  terms  f and  g and  the  actual  kinetics  will 
necessarily  be  obscured.  In  any  case,  the  present  analysis  does  shed 
light  on  the  role  which  activation  and  inhibition  mechanisms,  with  distinct 
migration  rates,  have  in  inducing  stability  of  nonuniform  structures. 
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We  deal  with  a two-component  system 
(20a)  u = u + f(u,  v) 

t XX 

(20b)  v = kv  + g(u,  v)  . 

t xx  ’ 

It  has  been  nondimensionahzed  so  that  the  diffusion  coefficient  of  u 
is  1;  we  let  k be  that  of  v. 

Definition:  u is  an  activator  if  f^  >0,  g^  > 0,  for  the  ranges  of  u 

and  v at  hand,  v is  an  inhibitor  if  f < 0.  q <0  for  this  ranqe. 

v v 

Thus,  increasing  the  amount  of  the  activator  increases  the  rate 
of  production  of  u and  v,  whereas  increasing  v has  the  opposite 
effect. 

First,  we  deal  with  "single  peak  distributions".  Let  «»(x),  ji(x) 

be  a pair  of  c"1  functions  which,  as  in  Fig.  5,  are  even  in  x,  approach 

limits  as  lx  I - «,  satisfy 

«>'(x)  * 0 for  x * 0,  <p"(0 ) * 0,  ip  = h(ijj) 

for  some  function  h with  h'(v)  >0  in  the  closure  of  the  range  of  4>> 

and  such  that  ^ . — : approaches  a nonzero  finite  limit  as  I x | — 

<p(x)  - <p(°°) 

Theorem  1:  Let  (<*>,^)  satisfy  the  above  assumptions.  Let  k > 1.  Then 

there  exist  functions  f(u,v),  g(u,v),  and  p(v)  such  that  f > 0, 

g > 0,  f < 0,  g < 0,  p > 0,  and  («/;,  ip)  is  a C°-stable  stationary 
u v v 

solution  of 

(21a)  ut  = uxx  + f(u>  v) 

(21b)  v = k(p(v))  (p(v)v  ) + g(u,  v)  . 

1 x x Q* 
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X 

Fig.  5.  Typical  single-peak  distribution. 

Moreover  if  <p  = c.^  4 c.,  for  some  constants  c with  c,  > 0,  then 

12  l 1 

p s 1,  so  that  our  system  is  of  the  form  (20). 

Remark:  In  fact,  our  construction  will  yield  many  solutions  of  the  inverse 
problem. 

Sketch  of  proof:  First  we  consider  the  case  s «^».  Under  the  assumptions 
given,  there  is  a C function  F(u)  such  that  <p"  = -F(«>).  The 
simplest  system  for  which  {<?,  4*)  is  a pattern  will  be  of  the  form 
(22a)  u = u 4 F(u)  4 <r(u  - v) 

t XX 

(22b)  v = kv  4 kF(v)  4 kn(u  - v)  , 

t XX 

where  <t  is  a sufficiently  large  real  number.  Other  examples  will  be 
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shuwn  later.  It  is  clearly  ol  activator-inhibitor  type  provided  is 
large  enough. 

It  is  seen  that  (v,  4>)  is  a solution  of  (22).  Wo  need  to  show 
it  is  stable.  Let  e = ~ < 1,  and  let  I.  be  the  linear  operator 
Lu  u'  + F'(v>(x))u.  Let  S be  the  linearization  of  the  right  hand  side 
of  (22)  about  (<p,  4*): 


Lu  + o (u  - v) 
k(Lv  + a(u  - v)) 

as  an  operator  from  C X C into  C x C . For  any  given  scalar  K, 
let  P (L)  be  the  operator 


PV(L)  £ t\z  f \[  (1  - e)n  - (1  + e)L | + LZ  , 


acting  on  functions  in  C (R).  If  V is  such  that  P (L)  is  boundedly 
invertible  on  C (R),  then  it  turns  out  that  (S  - \)  is  boundedly 
invertible  on  C°(R)  X C°(R),  so  that  \^Z(S)  (spectrum  of  S).  This 


is  proved  by  actually  writing  down  a solution  of  (S  - V 


terms  of  P (L)”1  and  L.  Hence  Z(S)  C {\  : z(P  (L))  3 0 },  so  that 
information  about  the  spectrum  of  S can  be  found  from  knowledge  of 


in 


Z(Px(L))  = PX(Z(L)). 

The  following  facts  about  the  function  P (p.)  can  be  verified 
directly:  For  fixed  p,  it  is  a quadratic  polynomial  in  Given  any 
M,  its  two  roots  will  be  real  and  negative  for  any  p < M and  p * 0, 
as  long  as  0 < e<  1 and  (1  - e)o  is  large  enough  (depending  on  M) 


-30- 


(Incidentally,  this  is  where  the  condition  k > 1 is  seen  to  be  necessary.) 

If  p - 0,  however,  one  of  the  roots  will  be  zero  and  the  other  negative. 

Now  the  spectrum  Z(L)  is,  in  fact,  real,  bounded  from  above  by  some 
number  M,  and  contains  0 as  a simple  isolated  eigenvalue.  This 
follows  from  our  assumptions  on  <p,  and  from  known  facts  about  the 
spectra  of  second  order  operators.  Thus  as  p ranges  over  Z(L),  the 
set  of  roots  i of  P (p)  will  range  over  a set  of  nonpositive  real 

A 

numbers  which  contains  the  origin  once,  and  is  otherwise  bounded  away 
from  zero.  In  view  of  the  fact  (shown  in  the  previous  paragraph)  that 
Z(S)  A,  we  know  that  Z(S)  has  the  same  property. 

But  with  this  knowledge  about  the  spectrum  of  S,  we  can  apply 

, , 0 

Sattinger's  stability  tneorem  | <14  | to  deduce  the  C -stability  of  the 
solution  u = <p,  v = <p.  (Incidentally,  I believe  this  to  be  one  of  the 
first  examples  in  which  Sattinger's  theory  has  been  applied  to  specific 
systems  of  order  greater  than  one.  ) 

Now  consider  the  case  when  v hfvj,)  £ if/.  We  write  the  system  (22) 
with  v replaced  by  the  symbol  w,  then  transform  the  second  dependent 
variable  according  to  the  equation  w h(v),  to  obtain  a system  of 
the  form  (21),  with  p(v)  h'(v).  Clearly  if  h'  = const,  we  may  take 
p 1.  This  completes  the  proof. 

Next,  we  consider  distributions  with  many  peaks,  in  fact  those 
which  are  periodic  in  x.  Let  «r( x)  and  4i(x)  be  periodic  C*’  functions 
such  that  <p"  = -F(<p)  for  some  function  F,  and  <p  = c^4>  + c^.  > 0. 
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Theorem  Let  «.  and  o satisfy  these  assumptions.  Let  k > 1. 

Then  there  exist  functions  f,  g satisfying  the  inequalities  stated  in 

Theorem  1,  such  that  (<.  , 4>)  is  a solution  of  (20).  Moreover,  the 

spectrum  of  S (the  linearization  of  the  right  hand  side  of  (20)  about 

u - e v - 4.)  is  contained  on  the  nonpositive  real  axis. 

Remark:  The  fact  tnat  none  of  the  spectrum  of  S is  in  the  unstable  half- 

plane  is  an  indication  of  the  stabilizing  influence  of  the  activator- 

inhibitor  mechanism,  although  it  is  not  proved  that  (<*>,  4>)  is  stable  in 
0 

the  C sense. 

The  proof  of  Theorem  2 proceeds  as  did  the  other  theorem,  with 
the  difference  that  it  is  no  longer  true  that  0 will  be  an  isolated  point 
of  the  spectrum  of  L. 

Qtner  solutions  of  the  inverse  problem:  Our  solutions  of  the  inverse 
problem  (in  the  case  <p  = 4j)  have  been  of  the  form  (22).  Equally  valid 
solutions  can  be  obtained  by  replacing  the  constant  a in  the  second 
equation  (22b)  by  a different  large  constant  r in  a certain  range 
depending  on  a,  and  by  adjoining  extra  terms  (u-v)  f(u,  v)  and 
(u  - v)^g(u,  v)  to  the  right  sides  of  (22a)  and  (22b),  respectively, 
where  f and  g are  arbitrary.  The  given  distributions  will  of  course 
still  be  solutions  of  the  revised  equations,  and  the  stability  analysis 
is  unchanged.  Moreover,  the  inequalities  given  in  the  statement  of  the 
theorem  hold  true  for  u near  <£,  v near  4/,  but  are  no  longer 
necessarily  true  for  arbitrary  u and  v. 
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beX/d 

Example:  The  function  t(x)  Tm~>  1S  an  even  peak-like 

(1  + e 7 ) 

distribution  with  < (0)  = 3/2,  </'(00)  = 0,  and  with  peak-width  of  the 

-2  2 

order  d.  It  can  be  checked  that  w satisfies  the  equation  <p"  + d (<p  - <p 
Setting  o(x)  - a<?(x),  a > 0,  we  see  that  (<?,  will  be  an  activator- 
inhibitor  pattern  for  the  system 

u = u + d 2(u2  - u)  + cr(u  - av)  + (u  - av)  f(u,  v)  , 
t xx 

v - kv  + kd  (av  - v)  + cx(u  - av)  + (u  - av)  g(u,  v)  , 
t xx 

provided  that  k > 1,  and  <r  is  sufficiently  large.  The  functions  f 
and  g are  arbitrary. 


A 
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III.  B.  Z.  Tlateaus 


The  previous  section  is  relevant  when  one  wishes  to  model  the 
concentration  o(  reacting  substances  at  isolated  locations.  The  present 
section,  on  the  other  hand,  is  concerned  with  modelling  phenomena  of 
difrerentiation,  characterized  by  the  development  of  sharp  boundaries 
between  regions  within  which  the  concentrations  are  relatively  uniform. 
We  call  these  sharp  boundaries  "transition  layers"  (see  [9]). 

The  methods  of  singular  perturbations  may  be  used  to  construct 
patterns  with  sharp  transition  layers.  This  was  explained  in  [ 10  J for 
boundary  value  problems  on  finite  domains  in  higher  dimensional  space, 
and  an  example  was  given  in  [9]  of  transition  layer  solutions  in  one 
independent  variable,  though  no  discussion  of  stability  was  given  in  the 
latter  paper. 

Here  we  review  the  construction  process  for  periodic  plateau-type 
patterns  on  the  whole  line.  In  this  case,  contrary  to  that  studied  in  [ 10 
there  will  be  no  boundary  effect  on  the  appearance  of  our  dissipative 
structures,  because  there  is  no  boundary.  The  following  will  be  formal 
and  nonrigorous . 


As  in  [ 10  |,  we  take  as  model  system  the  following: 

(23a)  u = e^u  f f(u,  v) 

t xx  ’ 


Fig.  6.  Regions  01  positivity  and  negativity  for  the  source 
function  f(u,v). 


Our  object  will  be  to  construct  periodic  stationary  solutions  having 

the  appearance  shown  in  Fig.  7.  As  seen  there,  the  upward  transitions 

occur  roughly  at  locations  a with  a = + np,  where  p is  the 

n n 0 

wave-length;  and  the  downward  transitions  are  at  (3^  = |3Q  + np.  On 

the  relatively  "Hat''  portions  between  the  transitions,  (2  3a)  with  u - e - 0 
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should  be  satisfied  approximately . That  is,  f(u,  v)  - 0,  so  the  phase 

plane  image  should  lie  near  the  S-shaped  curve  in  f ig.  6.  Referring 

to  that  figure,  where  the  two  ascending  branches  u - h (v)  are  attractors 

tor  the  c responding  kinetic  equauons,  we  specify  that  the  lower  flat 

regions  d < x < a be  mapped  onto  the  branch  u = h (v),  and  the 
n n+1  - ’ 

upper  portions  a < x < a be  mapped  onto  u = h (v). 

n n + 


Fig.  7.  Plateau-type  pattern  for  the  u-component. 

We  are  seeking  solutions  for  which  v does  not  change  abruptly 
through  the  transition  layer,  although  u does.  Therefore  to  lowest 
order,  there  is  a well-defined  value  of  v at  each  transition  point.  It 
is  explained  in  [ 10  | that  the  only  value  of  v which  can  support  a 
transition  layer  is  a value  v*  for  which 
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J(V*)  = 0 


f24) 

h+(v) 

where  J(v)  s f f(u,  v)du.  Let  us  assume  there  is  only  one  value 
h_(v) 

of  v satisfying  (24).  The  explanation  of  (24),  as  well  as  the  construction 
of  the  solution  inside  the  layer  itself,  is  along  the  following  lines. 

Suppose  the  transition  occurs  at  x = 0,  where  v = v . We  stretch 
variables  by  setting  | = x/e,  then  (23a)  with  u = 0 becomes,  to 
lowest  order  in  the  new  variables, 


(25) 


~ + f(u,  v ) - 0 . 

d£ 


The  methods  of  matched  asymptotics  dictate  that  (25)  should  be  solved 
for  £,  e (-or,*)  with  boundary  conditions  (if  it  is  an  upward  transition) 
u(-ao)  = h ( v ) , u(  + °°)  = h+(vQ).  Assuming  a solution  u(|)  exists, 
multiply  (23)  by  and  integrate  the  result  from  -°°  to  + <®  to  obtain 


0 = f 


oo  r 2 

I 4-  f l II  w , 

O'  d| 


1 _d  1 

'1  du) 

2 d|  1 

1 d£/ 

du 


h4(V0) 


de  = / 


h-(vo} 


f(u,  vQ)du  = J(v0)  ; 


hence  vy  vi;:.  It  turns  out  that  this  is  a sufficient,  as  well  as  a 
necessary,  condition  for  the  existence  of  a solution  of  (25)  with  the 
given  boundary  conditions. 

At  this  point,  we  have  obtained  the  following  information.  On  the 
broad  tops  and  bases  of  the  plateau  configurations  in  Fig.  7,  u and  v 
are  related  by  u = h±(v).  At  the  transition  points  between  them,  v = v*. 
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We  now  proceed  to  a determination  of  v(x).  In  (2  3b),  we  replace  the 
variable  u by  h^v),  depending  on  whether  we  are  on  the  top  or  the 
base  of  a plateau.  We  thus  ignore  the  fine  structure  c.'  the  transition, 


replacing  it  by  an  actual  discontinuity.  Then  (23b)  '.with  v^  - 0)  becomes 


(26) 


v + G(v,  x)  = 0 

XX 


where 


G(v,  x)  = 


g_(v)  = g(h_(v),  v)  for  Pn  < x < an+1  > 
g+(v)  s g(h+(v),  v)  for  a < x < P • 


Tiiis  equation,  together  with  the  condition  v = v*  for  x = or  (3^, 

can  be  expected  typically  either  to  have  no  solution,  or  a one-parameter 

family  of  periodic  solutions  (translations  in  x would  give  a second 

parameter,  which  we  ignore).  For  example,  suppose  g (v)  < 0.  Then 

v >0  for  x ( 1(5  ,a  ,),  which  implies  v < v*  for  x e (R  , a ,) 
xx  'Kn’  n+1  vrn  n +1 

and  v (0  ) < 0,  v (a  ) > 0.  If  g ,(v)  <0  as  well,  then  the  same 
x n x n + 

derivative  conditions  must  hold  at  the  endpoints  of  the  intervals  (a  , i3  ). 

n n 

But  these  two  sets  of  derivative  conditions  contradict  one  another,  so 
there  can  be  no  solution. 

For  a pattern  to  exist,  it  is  therefore  necessary  that  g±(v)  not 
be  of  one  sign.  As  a simple  example,  let  us  assume  that  g±(v)  are 


constants  with  g < 0 < g + . Then  from  (26),  we  have  v"  = -g  ^ 

x t (a  , t3  );  hence  setting  p , = (3  - a , we  get 

n n + n n 


for 


v = V.  * - gt(*  - «n)(Hn  - X),  X < (V  dn)i  v;<un> 


g p 
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Thus  in  this  simple  example,  we  apparently  obtain  a plateau-like  solution 
for  any  spacing  p of  the  o's  and  (Vs  satisfying  (27).  However, 
this  is  only  partly  true:  the  wavelength  p cannot  be  too  large,  for 
then  the  range  of  the  function  v might  extend  beyond  the  domains  of 
one  or  both  of  the  functions  h . 

The  foregoing  describes  a method  of  obtaining  formally  approximate 
solutions  for  small  e.  The  procedure  could  be  continued  to  yield  higher 
order  approximations;  the  constructions  in  [10)  and  [11]  are  relevant  here. 

Now  supposing  a plateau  solution  (which  we  denote  by  (u^(x),  v^(x))) 
has  been  constructed,  how  can  we  infer  its  stability  or  instability?  The 
following  heuristic  argument  suggests  a workable  criterion. 

Suppose  we  subject  the  given  stationary  solution  to  a perturbation 
which  is  small  in  the  C norm.  Such  a perturbation,  of  course,  does 


not  erase  or  relocate  the  transition  layers.  For  the  moment,  let  us 
assume  that  the  latter  do  not,  in  fact,  move  during  the  subsequent  evolu- 
tion of  the  perturbed  system.  Then  between  the  layers,  the  image  points 
(u,v)  will  again  be  drawn  to  the  attractors  h±(v)  (by  neglecting  e 
in  (23a)),  and  at  the  same  time,  the  evolution  of  v will  be  governed 
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by  t23b).  Eventually,  the  term  g in  this  latter  equation  may  be 

replaced  by  G(v,x)  (see  (26)): 

i28)  v = v + G(v,  x)  . 

t xx  ’ 

The  solution  before  perturbation  was  such  that  v satisfied  (28)  and 
was  time-independent,  of  course.  So  we  must  now  investigate  the 
stability  of  solutions  of  (28).  If  we  impose  the  condition 

(29)  ^9±(v)<0 

(this  being  the  v-denvative  along  the  ascending  branches  in  Fig.  6), 
then  the  solution  vQ(x)  will  be  stable.  In  fact,  for  small  positive 
constants  6,  vQ(x)  ± 6 will  be  super-  and  sub-solutions  of  (28),  so 
again  as  in  [ 1),  solutions  of  (28)  with  initial  data  caught  between 
vQ  - 6 and  vQ  + 6 will  converge  uniformly  to  v^ . 

This  argument  asserts  stability  under  condition  (29)  if  we  assume 
the  transition  layers  are  immobile.  They  are  not,  in  fact,  immobile, 
but  their  movement  is  slow  and  can  be  analyzed  by  the  techniques  in  | 11 1 . 
Immediately  following  the  initial  perturbation,  the  value  of  v at  any 
given  transition  layer  will  no  longer  necessarily  satisfy  J(v)  = 0 (as 
in  ( 24)).  What  this  means  is  that  the  layer  will  develop  into  a sharp 

X — GCt 

wave  front  which,  locally,  will  be  of  the  form  u = U(  - ,x),  where 

c depends  on  v at  the  position  of  the  front,  so  will  itself  change  with 
time.  This  approximate  representation  of  u follows  the  argument  in 
| 10  |,  | 11 1,  to  which  the  reader  is  referred  for  more  details.  As  stated 
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above,  the  velocity  ec  will  be  small,  so  that  until  a time  of  the  order 
1/e  has  elapsed,  the  above  assumption  that  the  layer  is  immobile  should 
be  valid. 

But  after  this  time,  shifts  in  the  transition  layer's  position  due  to 
the  slow  motion  of  the  front  may  have  accrued. 

The  effect  of  changes  in  the  layer's  position  should  therefore  be 
determined.  It  turns  out  that  the  nature  of  this  effect  depends  on  the 
signs  of  the  two  functions  g±(v).  In  the  following,  we  call  the  collection 
of  intervals  on  which  u ~ h+(v)  the  domain  of  the  (+)  state,  and 
similarly  for  the  (-)  state.  Suppose,  for  example,  that  g (v)  < 0 < g (v). 

It  will  then  also  be  true  that  g(u,  v)  < 0 when  the  point  (u,  v)  is  in  a 
neighborhood  of  the  left  ascending  branch  in  Fig.  6,  and  g > 0 near  the 
right  ascending  branch.  We  focus  attention  on  the  effect  of  shifts  in 
the  position  of  a single  transition  layer.  Suppose,  initially,  that  we  are 
in  a stationary  state,  and  that  the  points  (u,  v)  lie  in  the  two  neighborhoods 
described  above  on  opposite  sides  of  the  layer.  Now  suppose  the  position 
of  the  layer  is  perturbed  in  the  direction  which  will  increase  the  "domain" 
of  the  state  (+)  (that  is  to  say,  it  is  moved  away  from  the  ( +)  state 
toward  the  (-)  state.  Since  gf  > g , the  immediate  effect  of  this  is 
to  increase  the  source  term  g(u,  v)  in  (23b),  in  the  interval  between  the 
old  and  the  new  positions.  This  will  cause  the  function  v^x,  t),  hence 
v,  to  increase.  At  the  same  time,  the  value  of  v at  the  layer  has 
also  changed  by  virtue  of  the  fact  that  the  initial  stationary  function  v(x) 
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was  not  constant,  whereas  the  layer's  position  was  altered.  This  second 


effect  may,  for  a period  of  time,  counteract  the  first  effect  (the  increase 
of  v due  to  increasing  the  source  terms);  but  eventually  the  first  will 
dominate,  because  the  subsequent  motion  of  the  layer  is  slow,  and  v( 
will  remain  positive  until  the  position  has  moved  a significant  amount. 

So  eventually  the  value  of  v at  the  layer  will  surpass  its  initial  value. 

Let  us  now  assume  that  J(v)  is  an  increasing  function  of  v 
(apparently  the  case  for  the  function  f depicted  in  Fig.  6).  Then  such 
an  increase  in  the  value  of  v at  the  front  also  raises  J,  which  in  turn 
results  in  the  (+)  state  becoming  more  dominant,  in  the  terminology  of 
| 10  | and  [ 11 ) . This  means  that  the  (slow)  motion  of  the  front  will  proceed 
in  a direction  which  further  increases  the  domain  of  the  (+)  state.  The 
perturbed  layer  will  therefore  not  tend  to  return  to  its  original  location, 
and  we  are  in  an  unstable  situation. 

On  the  other  hand,  if  we  assume  the  opposite  inequalities 
(30)  g+(v)  < 0 < g_(v)  , 

then  the  above  analysis  indicates  that  a perturbed  position  of  the  layer 
will  eventually  move  back,  toward  its  original  location.  Essentially  the 
same  analysis  holds  if  the  perturbation  is  in  the  opposite  direction. 

On  the  basis  of  the  above  arguments,  which  admittedly  are  far  from 
being  rigorous,  I conjecture  that  the  constructed  periodic  plateau-type 
stationary  solution  will  be  stable  if  (29)  and  (30)  hold,  and  if  J is  an 
increasing  function  of  v. 


In  [ 10  | , an  argument  for  stability  or  "realizability"  was  given,  in 


some  respects  similar  to  the  above.  Yet  in  other  respects,  it  was 
different,  because  the  effects  of  boundary  conditions  in  a bounded  domain 
were  essentially  involved. 


I 
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IV.  Discussion 


This  paper  is  devoted  entirely  to  patterns  on  the  entire  real  line, 
and  it  may  well  be  objected  that  real  reacting  and  diffusing  systems  are 
always  finite.  Suppose  the  real  system  is  one-dimensional,  and  its 
domain  is  a finite  interval  I.  If  the  length  of  I is  large  compared 
with  the  characteristic  length  of  the  patterns  formed,  then  the  infinite-line 
approximation  should  be  a good  one,  in  the  interior  of  I.  In  fact,  the 
patterns  studied  here  are  generated  by  the  processes  of  reaction  and 
diffusion  only,  and  so  should  exist  in  interior  regions  independently  of 
any  boundary  effects.  Of  course  there  also  exist  patterns  which  are 
generated  or  strongly  conditioned  by  the  presence  of  boundaries  and  by 
the  transport  mechanisms  occurring  there.  This  type  of  pattern  is  not 
subject  of  the  paper.  In  short,  our  restriction  to  the  infinite  line  was 
made  for  the  purpose  of  exploring  pattern  phenomena  generated  solely  by 
reaction  and  diffusion,  without  having  to  worry  about  boundary  effects. 

Nevertheless,  some  things  can  be  said  immediately  about  certain 
patterns  on  finite  domains,  analogous  to  those  discussed  here. 

Consider  first  the  scalar  case,  n 1.  This  is  a matter  of  solving 
(Z)  on  an  interval,  with  given  boundary  conditions.  We  restrict  attention 
to  homogeneous  Neumann  (no-flux)  or  to  Dirichlet  (u  = 0)  boundary 
conditions.  Again,  all  solutions  can  be  constructed,  and  their  stability 
analyzed  by  techniques  similar  to  those  described.  It  turns  out  that  under 
Neumann  conditions,  the  only  stable  solutions  are  certain  constant  ones, 
and  under  Dirichlet  conditions,  the  only  stable  ones  are  those  which  have 
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a single  internal  extremum  (except  in  the  case  when  u 0 is  stable). 

See  also  the  treatment  in  ( =>]. 

Consider  next  the  question  of  small  amplitude  patterns,  as  treated 
in  III.  A.  Set  on  a finite  interval  with  homogeneous  boundary  conditions, 
this  is  a standard  bifurcation  problem,  and  the  results  of  [ 6 | , ( 7 j , [2  1] 
apply.  One  commonly  finds  stable  super-critical  bifurcation  branches 
of  nonconstant  solutions  which,  for  small  amplitudes,  approximate 
periodic  functions  with  many  wavelengths  fitted  into  the  given  interval 
([2],  1-4]).  These  will  also  approximate  some  of  the  periodic  solutions 
found  in  III.  A. 

The  treatment  of  peaks  — single  or  multiple — in  III.  B.  1.  may  be 
extended  with  little  change  to  the  case  of  a finite  interval  with  (say) 
Neumann  boundary  conditions  imposed.  (Of  course,  one  must  begin  with 
distributions  v1  and  ^ which,  themselves,  satisfy  these  boundary 
conditions).  The  same  is  true  of  the  plateaus  in  III.  B.  2.  See  especially 
[9]  and  [10]. 

A second  objection  may  be  that  real  reacting  and  diffusing  systems  are 

distributed  in  space  or  on  a surface,  rather  than  on  a line,  as  we  have 

here.  It  is  probable  that  the  extension  of  the  present  results  to  problems 
2 3 

on  R or  R will  entail  many  complications  and  difficulties.  In  the 
case  n 1,  for  example,  one  can  no  longer  immediately  write  down 
all  stationary  solutions  of  Au  + f(u)  = 0;  and  though  in  principle  the 
technique  cJ.Arnnson  and  Weinberger  works  in  space,  its  application  may 
be  difficult. 
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As  for  small  amplitude  patterns  in  the  case  n > 1,  this  problem  is 
frauyht  with  most  of  the  difficulties  associated  with  the  Benard  problem 
in  fluid  mechanics.  One  can  find  solutions  which  are  doubly  periodic 
with  various  kinds  of  symmetries,  but  the  stability  question  is  very 
difficult,  unless  one  allows  only  perturbations  within  certain  classes 
restricted  by  artificial  symmetry  conditions.  Of  course,  if  one  wants  to 
find  patterns  for  which  boundary  effects  are  important,  one  has  a boundary 
value  problem  in  a bounded  domain.  Typically,  this  reduces  to  a standard 
bifurcation  problem,  and  the  identification  of  the  bifurcation  points, 
together  with  the  determination  of  the  approximate  shape  of  the  small 
amplitude  patterns  which  arise,  involves  solving  a linear  eigenvalue 
problem . 

There  are  real  complications  involved  in  extending  the  treatment  in 
III.  B.  1 to  higher  dimensions . On  the  other  hand,  plateau  configurations 
(as  in  III.  B.  2)  can  be  constructed,  and  was  done  in  ( 10]  (see  also  [ 12]) 

A final  comment  should  be  made  about  the  terminology  "pattern 
formation".  The  same  term  is  also  commonly  used  in  developmental 
biology,  where  it  refers  to  structural  organization  in  developing  organisms 
The  possible  relevance  of  reaction-diffusion  patterns,  such  as  those  in 
the  present  paper,  to  patterns  in  this  other  sense,  will  suggest  itself. 

No  claim  is  being  made  here  that  reaction-diffusion  patterns  are  anywhere 
near  to  being  adequate  models  for  the  biological  patterns.  In  fact,  the 


latter  are  undoubtedly  highly  complicated  phenomena  involving  many  other 
types  of  processes  besides  those  suggested  here.  At  the  same  time, 
any  serious  student  of  theoretical  developmental  biology  should  be 
aware  of  the  rich  repertoire  of  solutions,  including  stable  patterned  ones, 
exhibited  by  rather  simple  models  based  on  reaction  and  diffusion  alone. 
To  this  extent,  I hope  the  results  given  in  this  paper  may  be  useful  to 
the  biologist. 
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